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ABSTRACT 

Multiwavelength data indicate that the X-ray emitting plasma in the cores of galaxy clusters is not 
cooling catastrophically. To large extent, cooling is offset by heating due to active galactic nuclei 
(AGN) via jets. The cool-core clusters, with cooler/denser plasmas, show multiphase gas and signs of 
some cooling in their cores. These observations suggest that the cool core is locally thermally unstable 
while maintaining global thermal equilibrium. Using high-resolution, three-dimensional simulations we 
study the formation of multiphase gas in cluster cores heated by highly-collimated bipolar AGN jets. 
Our key conclusion is that spatially extended multiphase filaments form only when the instantaneous 
ratio of the thermal instability and free-fall timescales (iTi/iff) falls below a critical threshold of « 10. 
When this happens, dense cold gas decouples from the hot ICM phase and generates inhomogeneous 
and spatially extended Ha filaments. These cold gas clumps and filaments 'rain' down onto the central 
regions of the core, forming a cold rotating torus and in part feeding the supermassive black hole. 
Consequently, the self-regulated feedback enhances AGN heating and the core returns to a higher 
entropy level with frri/iff > 10. Eventually the core reaches quasi-stable global thermal equilibrium, 
and cold filaments condense out of the hot ICM whenever in/iff % 10. This occurs despite the fact 
that the energy from AGN jets is supplied to the core in a highly anisotropic fashion. The effective 
spatial redistribution of heat is enabled in part by the turbulent motions in the wake of freely-falling 
cold filaments. Increased AGN activity can locally reverse the cold gas flow, launching cold filamentary 
gas away from the cluster center. Our criterion for the condensation of spatially extended cold gas is 
in agreement with observations and previous idealized simulations. 

Subject headings: cooling flows - galaxies: active - galaxies: clusters: general - galaxies: jets - 
intergalactic medium - methods: numerical 



1. INTRODUCTION 

Unopposed radiative cooling of the intracluster 
medium (ICM) would lead to the cooling catastrophe 
- an inflow of gas toward the center, due to the loss 
of pressure, a t the rate of several hundred M yr _1 (see 
Fabian|1994 for a review) . The main consequences would 



be strongly peaked surface brightness maps, central tem- 
peratures falling below T v - ir /3, and a huge accumulation 
of cold gas in the cluster core, well over 10 12 M Q . The 



(e.g., Peterson et al. 


2001 


2003 Tamura et al. 12003 Pe- 


terson & Pabian^OO 


i| points out that either the cooling 



of heating. 

The most popular and promising heating mechanism 
is the active galactic nuclei (AGN) feedback. Radio and 
X-ray data provide overwhelming evidence for the in- 
teraction of jets/outflows from t he supermassive black 
holes with the surrounding ICM flBohringer et al.l|1993| 
Blant on et al.|[200Tj |Finoguenov fc J ones 2001; Jones 



et al. 



2007 



2002; IRafferty et al. 2006; JMcNamara fc Nulsen 



Finoguenov et al.||2008| |Rafferty et al.||2008| ). The 



power supplied by the AGN is in the range of 10 4ci — 10 46 
erg s _1 , adequate to globally offset radiative cooling 

6 Work performed while on leave at the Depertment of Astron- 
omy in the University of Michigan, Ann Arbor. 



2003 



2007 



losses. Further support for the outflow model comes from 
blueshifted absorption l ines in the core of several g alax- 
ies: in the optical band (Ncsvadba et al. 2008, 2011|), UV 
and X -ray bands (Tombesi et al. 2010; Crensha w et al 
for a review), and at 21-cm ( pYlorganti et al.potl" 
. These lines originate at kpc distances and corre 
spond to velocities of up to several 10 4 km s _1 and mass 
outflow rates around a M Q yr _1 or more. The geometry 
is still unclear, but radio images are suggestive of con- 
ical outflows with narrow opening angles and possibly 
consistent with precession. 

While the global energy balance may be achieved with 
the he lp of the AGN feedback and/or other mecha- 
nisms (|Ru szkowski & Bcgclman 2002 ; Brighcnti & Math 
lews] 120031 lUuo et al. 1 12008| |Conroy fc Ustrrker||200~ 
Ruszkows ki fc Oh||201ip , the issues regarding the cou 
pling of the mechanical AGN energy to the thermal ICM 
energy, along with the spatial distribution of heating, are 
not settled. The coupling may occur, for examp le, via 
turbulent mixing of AGN bubbles with the ICM (Scan- 
napieco & Briiggcn 2008), viscou s dissipation of weak 



shocks (Ruszkowski et al. 2004a b|, and cosmi c ray heat 
ing (|Guo fc Oh||2008| IRuszkowski et al.||2008| |Sharma et 



|al.|20O9) . Other mechanisms may also help heat the ICM 
and redistribute th e energy throughout the cool core: 
dynamic al friction (lELZant et al. 2004) , thermal con- 
duction (JNarayan fc Medvedcv 2001; Bogdanovic et al. 
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Gaspari, Ruszkowski & Sharma 



2009] |Ruszkowski fe Oh|2010l [201ll [Parrish et al.1[20lot, 

or slo shing motions (Morso ny et al.p OlO; ZuHone et all. 
20101. 



Tot most microscopic heating mechanisms the ICM 
is expected to be locally thermally unstable. In other 
words, even if the ICM maintains global thermal bal- 
ance without a cooling flow, some of the gas can ex- 
perience run awa y cooling. Indeed, t he observations by 



Edge| ( |2001| ) and |Salome et al.| ( |2003| ) revealed that cool- 
ing flows are not totally stopped: the mass of cold gas, 
mainly dominated by molecular H2 associated with CO 
emission, reaches 10 9 — 3 x 10 11 M Q in the majority of 
cool-core clusters. Cooling may also be inferred from 



star formation (e.g., O'Dea 2008, typically 1-10 M, 
1 ' More r ecently, a wide survey conducted by Mc 



yr 



Donald et al] (2010, 20 11a|b[ ) utilized Ha emission as a 
robust tracer of the cold phase. In more than 60% of 
clusters and groups, extended and/or nuclear cold gas 
was detected. The morphology of this phase is usually 
clumpy or filamentary and e longated toward the center 



(e.g., in the Perseus cluster; Salome et al. 2006 Fabian 



et al. 2008). 



aid et al. 



Cavagnolo et al . 



h 



2008|), Voit et al. (|2008|) and McDon- 
tne" 



ave shown the existence of a tight 
correlation between the amount of cold gas and the ICM 
entropy or the cooling time. Only clusters with central 
entropies below 30 keV cm 2 show evidence for the ex- 
tended cold gas. This criterion can be expressed in terms 
of the ratio of the cooling time (i C ooQ and the free-fall 
time (is). 



Analytical e stima tes and Cartesian simulations by Mc- 



Court et al. (2011 



Mil) pointed out that a spatially 
extended and inhomogeneous cold phase can form when- 
ever t cooi/tg < 1. Previous studies (e.g., Balbus & Soker 
19891 did not predict the presence of multiphase gas 
because they neglected a spatially distributed heating, 
which provides global thermal stability. Fig. 12 in Mil, 
based on observational data, sugg ests that the c r itical 
threshold is actually higher, ss 10. Sharma et al. (2011 
- Sll) extended the work of Mil, considering spherical 
feedback models. They found that, due to the geomet- 
rical compression, the criterion is in fact less stringent: 
instability occurs when t CQO \/tg < 10. 

The results presented in Mil and Sll are very promis- 
ing, but the distributed heating models adopted there are 
idealized. In these models the heating is either added us- 
ing a prescribed spherically-symmetric functional form, 
or it is constructed to be equal to the shell-averaged ra- 
diative cooling loss at every radius. Given the simplicity 
of these prescriptions, it is necessary to consider more 
realistic heating models to reliably study the formation 
of the multiphase medium, along with reproducing sev- 
eral important features, such as shocks, cavities and tur- 
bulence. We perform high-resolution three-dimensional 
simulations which incorporate ICM heating via colli- 
mated AGN jets. This 
of our previous efforts 
Glla,b), which included 
focused on the dynamics of the hot ICM phase, suppress- 
ing the formation of the multiphase gas via cold phase 

2 The cooling time and the thermal instability timescale are 
similar for a microscopic heating mechanism which is independent 
of density, so we use txi and t C ool interchangeably. 
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2009 




2011a|b 
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dropout. 

Our more realistic AGN feedback simulations confirm 
that extended multiphase ICM forms only when the 
instantaneous ratio of the thermal instability and free- 
fall timescales falls below a critical threshold of w 10, 
irrespective of the initial condition. The cluster core 
attains rough thermal equilibrium after a few cooling 
times. In addition to occasional extended cold filaments, 
centrally concentrated cold gas exists at most times. 
Thus AGN feedback primarily occurs via the cold phase. 
Our simulations also show that sometimes the cold gas 
can be carried out by AGN jets. The evidence for such 
cold filaments dredged out from the central galaxy is 
seen in the Perseus cluster ( Fabian et al.||2003| . 



2. NUMERICAL METHODS 



We use the adaptiv e mesh refinement code FLASH 4.0 



(Fryxell et al. 2000), modified to simulate mechanical 
AGIN feedback m the cooling ICM. The numerical finite- 
volume method employed is the third-order accurate split 
PPM, suited to solve the Euler equations including the 
required heating and cooling terms: 



dp 
dt 



V • (pv) = ±V • (TO jet n z ), 



(1) 



dpv 

~dt 



-V-(pv® v) + VP = /9g±V-(TO jct n z ® v jct ), (2) 



dps 
~dt 



V • [(pe + P) v] = pv • g - n c niA(T, Z) 



P=( 7 -l)p 



± V ■ ( -mjotVjet n z ) , (3) 



(4) 



where p is the gas density, v the velocity, e the specific 
total energy (internal plus kinetic), P the pressure, and 
7 = 5/3 the adiabatic index. The temperature is com- 
puted from P and p using Eq. Q, with an atomic weight 
(i ~ 0.62, appropriate for a totally ionized plasma with 
25% He in mass. The externally imposed gravitational 
acceleration g accounts for the dark matter NFW halo 
( Navarro, Frenk fc White|1996 ) characterized by a virial 
mass and radius of 10 io M and 2.6 Mpc, respectively 
(concentration parameter ~ 6.6). Our implementations 
of the cooling and jet flux terms in Equations |l]) - ([3| 
deserve careful discussion and are presented in Sections 
2.1 and 2.2, respectively. 

We decided to neglect magnetic fields and anisotropic 
thermal conduction. Mil showed that the anisotropic 
conduction changes the morphology of the cold gas 
(more filamentary), but does not affect the criterion for 
the generation of the multiphase gas due to thermal 
instabilities. We defer the study of these effects to future 
work. 



2.1. Radiative Cooling 

The intracluster medium emits radiation mainly in the 
X-ray band due to Bremsstrahlung (T > 10 7 K) and due 
to line emission at lower temperatures (10 4 - 10 7 K). 
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Assuming collisional ionization equilibrium, the cooling 
source term in Eq. ^ can be modeled with a cooling 
function A(T,Z): we use an analytical fit to the Suther- 
land & Dopita ( 1993[ ) tabulated values for a fully ionized 



plasm a with metallicity Z = 0.3 Z (e.g., Tamura et al. 



2004 1 , setting the temperature floor at 10 4 K. The total 



emissivity of the gas is also proportional to the electron 
and ion number densities: C = n c riiA(T, Z). For simplic- 
ity we assume that molecular weights are constant, with 
HelH = 1.464. 

An important difference between the current imple- 
mentation of the cooling t erm and that adopted in our 
previous models (Glla,b; |Brighenti fc Mathews||2002| ) 
is the absence of the dropout term, which removes the 
cold gas from the domain. Here we are interested in the 
formation and growth of thermal instabilities, producing 
cold clumps and filaments in the ICM. 

The presence of the cooling term can significantly slow 
down the computation. Explicit schemes, such as Runge- 
Kutta, are usually simple and fast (per timestep), but 
impose severe restrictions on the timestep At, making it 
much smaller than the regular hydrodynamical timestep. 
Implicit schemes allow to use larger timesteps, but these 
methods are less precise. In a splitting framework, i.e., 
first solving the hydrodynamic PDE and then the cooling 
source term ODE (assuming fixed density), no implicit 
or explicit solver is actually needed for the ODE step. 

The new temperature T" +1 of a cell j at time n + 1 
can be computed from the old temperature T™ by exactly 
evaluating the following integral: 



dT = (7-1)// 
A(T') k B fi c ^im p 



Pj 



At 



(5) 



A practical way to compute the integral is to split the 
integration interval into two parts: one from T™ to T 1G f 

and the other one from T re f to T" +1 , using an arbitrary 

reference value (see Townsend 2009). This leads to 



ijm+1 _ y~ 



Y{Tf) + 



Tf A(T ref ) At 



T Iei a(2?) t cool 

where the cooling time and Y are defined as: 
e _ kBHemrrip Tf 

= a(tw) r T - dT' 

{ >~ T ref Jt A(T')' 



t. 



cool 



(6) 



(7) 



(8) 



where e = P — 1) is the internal energy density. The 
temperature T" +1 can be exactly computed if the recip- 
rocal of the cooling function is analytically integrable. 
Thus, we decompose the cooling function into a large 
number of piece- wise powerlaws in logarithmic tempera- 
ture intervals 



A(T) = A fc 



Th <T <Ti 



fc+i 



(9) 



and precompute at initialization the constants of inte- 
gration for each interval. The integral in Eq. (Tsl) can be 



computed analytically using this powerlaw decomposi- 
tion of the cooling function. The only effort in retrieving 
T" +1 comes from finding the correct interval in which 

the argument of Y^ 1 resides (Eq. (kSJ)). 

Aside from being exact (in the ODE step), this method 
is faster than either explicit or implicit methods. The er- 
ror from the piece-wise powerlaw approximation is neg- 
ligible, because the number of the temperature intervals 
can be made very large (matching the tabulated Suther- 
land & Dopita values), without large memory require- 
ments. 

In principle, our solver does not require a timestep lim- 
iter for stability. However, the coupling between the hy- 
drodynamics and source term integration is still approxi- 
mate. Thus we preferred to impose an upper limit on the 
cooling timestep (~ t coo i)- We note that this did not dra- 
matically increase the computational requirements. In 
order to further improve the coupling, we have also im- 
plemented a second order accurate temporal evolution in 
the splitting scheme, i.e., we alternate the order of the 
integration of the hydrodynamical (H) and source terms 
(S) as follows: H-S-S-H. 

We have compared the explicit calculation, setting 
a strong limiter (< 0.1 t coo \), with the 'exact' cooling 
run and we obtained equivalent results (at least in the 
regime of parameters relevant to this study and down 
to the scales comparable to few smallest resolution 
elements) . 



2.2. AGN Heating 

The heating module consists of two key elements: com- 
putation of mass accretion onto the central supermassive 
black hole, and the energy injection via collimated bipo- 
lar outflows. The quantity of accreted material is es- 
timated as the inflow through a spherical surface very 
near the center, r ~ 500 pc. Each cubical cell j con- 
taining a piece Aj of the spherical surface, contributes 

■Mace — Pj v rjAj to the accretion rate only if the ra- 
dial velocity v rj is negative (inflow). The accreted gas, 

M acCj At, is then removed from the corresponding zone. 

After computing the total mass accretion rate, we cal- 
culate the kinetic power of the bipolar outflows from 



Pjot = eM a , 



(10) 

where e is the (mechanical) efficiency of the feedback, 
which is a free parameter of our model. We note that 
the actual accreted gas should be ideally calculated a 
few gravitational radii from the black hole. This would 
require an extreme dynamical range, which is currently 
not possible due to hardware limitations. However, we 
can simply absorb any uncertainty in the actual accretion 
rate into our efficiency parameter e. If only a fraction / 
of the gas entering our sink region falls into the black 
hole, the real efficiency would need to be rescaled to e/f. 

The AGN mass, momentum, and energy are injected 
into the volume through the internal boundaries in the 
middle of the computational domain. The total cross- 
sectional area of the nozzle is A noz = (2 Ax) 2 , where 
Ax is the highest resolution in the domain. Both jets 
originate at the same interface along the positive and 
negative z-direction (n z is the z- vector). The jet flux 
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TABLE 1 

Model parameters 



Model 


efficiency e 


min(t coo i/t ff ) 


Notes 


pure CF 




7, 21 


no AGN heating 


r7-6em3 


6 X 1CT 3 


7 


along ±n z 


r7-6em3w 


6 x icr 3 


7 


random wobbling 


r7-6em3c 


6 x 10" 3 


7 


convergence (2 As) 


r7-lem2 


1CT 2 


7 


along ±n z 


r21-6em3 


6 x icr 3 


21 


along ±n z 


r21-lem2 


10" 2 


21 


along ±n z 



terms appear on the right hand side of Equations |T]), 
Q and This method of including AGN feedback 

appears to be slightly better than injecting quantities at 
zone centers, which may lead to artificial evacuation of 
the regions perpendicular to the jets. 

The integrated mass outflow rate, Mj ct = mjet^lnoz, 
is based on the conservation of momentum (Mj e tVj et ), 
leading to 



M jrt =2P iet /v 



jot: 



(11) 



where Vj e t is the subrelativistic outflow velocity fixed at 
5 x 10 4 km s _1 . The velocity of AGN outflows is still 
poorly constrained and observations sug gest values from 



several 10 3 to several 10 4 k m s 1 (e.g., 



2003 Tombcsi et al. 



velocity gives reasona 



Crenshaw et al. 



Our choice of the outflow 



2010J). 

)le values for the mass outflow rate 
of ~ Mq yr _1 . Notice that more massive (and slower) 
outflows tend to carve a more visible channel in the ICM, 
while lighter (and faster) jets release most of their energy 
in the central region, producing strong shocks and big 
bubbles similar to X-ray data (see Glla,b). 

We also implemented wobbling jets, with variable 
inclination (9) and azimuthal angle (0). In these 
computations the location of the nozzle remains fixed 
in space, but all three velocity components (e.g., 
v x = Vjet sm # cos0) vary with time. The injection is still 
"cylindrical" , in the sense that the velocity vectors at 
the nozzle remain parallel. We consider simple models 
where jet orientation varies randomly within a specified 
solid angle. These models are further discussed in 
Section 3.2.3. 



2.3. Initial Conditions 

We model our initial conditions on Abell 1795, which is 
a typical massive relaxed cool-core cluster characterized 
by the virial mass of ~ 10 15 M Q . In order to compute 
the initial hydrostatic gas density profile, we take 
the observed temperature fit (Fig. [II second column, 
black dashed line) and use the analytical formula for 
a NFW dark matter gravitational acceleration. The 
normalization of the density profile is chosen such that 
the gas fraction is s» 0.15 at the virial radius. The 
obtained initial conditions correspond to the initial 
minimum ratio of t coo \/tg ~ 7, where the free-fall time 
is given by 



2r 
9(r) 



1/2 



(12) 



we use outflow boundary conditions. We use static mesh 
refinement. The linear size of each shell doubles among 
adjacent levels. The highest refinement level is 10, and 
every block is divided into 8 3 zones. The finest zone of 
the box has a size of Ax ~ 300 pc. 

As noted by Sll, symmetric initial conditions do not 
produce thermal instabilities. Therefore, we perturb the 
initial density and temperature profiles to initiate the 
formation of cold filaments and blobs. Depending on the 
model parameters, such features can develop even in the 
early stages of evolution, before the AGN has become 
very powerful. We also note that the presence of the 
instabilities early on facilitates dispersal of the energy 
injected by the collimated AGN jets. 

The fluctuations leading to the thermal instability can 
be produced by the turbulence associated with cosmo- 
logical mergers. It is beyond the scope of this work to 
study such effects. It is thus sufficient to adopt a simple 
Gaussian random field with a white noise power spec- 
trum of gas density. Using dispersion relation and direct 
simulation, Mil and Sll demonstrated that thermal in- 
stability does not depend on the range of the driving 
modes in Fourier space. Moreover, in the late stage of 
the evolution, the density fluctuations are dominated by 
the AGN jets. 

We generate initial conditions using IDL. First, we ini- 
tialize a complex stochastic field (1024 3 ) in the k-space, 
with amplitude of each mode given by 



where k = (k 2 + k 2 



Wock~\ (13) 

+ k 2 ) 1 / 2 . The above prescription en- 
power spectrum, which is oc W 2 k 2 , 



sures that the '3D' 
remains constant, i.e., is described by a white noise spec- 
trum. Second, we combine the above amplitude wi th a 
Gaussian random field (e.g., |Ruszkowski et al.|[2007 ): 



[Re{W(k)), Im(W(k))] = [G(u x )W, G{u 2 )W] (14) 

where G is a function of a uniform random deviate u\ 
or U2 that returns Gaussian-distributed values. This 
method guarantees that the phases are random. We also 
impose cutoffs at fc max = 1024/4 and fc m i n = 3 (in units of 
2tt/L). The high-k cutoff is introduced in order to avoid 
spurious numerical effects. Third, we compute inverse 
Fourier transform of W(k) to obtain W(x, y, z). We nor- 
malize it by the mean of its absolute value (|W(x, y, z)\). 
Finally, we superpose the normalized fluctuations onto 
the initial density distribution p un p' 



p(x, y, z) = p unp (x, y,z)[l + £ 



W(x,y,z) 
■{\W(x,y,z)\) 



(15) 



where £ = 0.3 is the amplitude of the perturbations. 
The new temperature distribution is computed from the 
new density assuming isobaric fluctuations. 



We present 
efficiencies e, 
and wobbling jets) 

icool/ ts 



3. RESULTS 

simulation results for different AGN 
methods of energy injection (steady 



minimum 



Our computational domain is 1.27 Mpc on a side and 



and different initial values of the 
in the ICM (Tl-ratio; minimum 
value over all radii). Our previous results (Gila; Sll) 
demonstrated that a suitable choice of the mechanical 
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AGN efficiency may be in the range 10 3 
supported b y observational constraints 



10" 2 . This is 



e-g-, 



Merloni & 



Heinz|[2008 ) . Common observed values for the minimum 
initial Tl-ratio are in the range from 5 to 25 (Mil). 
Motivated by these results, we perform computations 
with initial Tl-ratio 7 or 21, and e = 6 x 1(T 3 or 1(T 2 
(see Table 1). 

3.1. Pure Cooling Flow 

We first analyze a reference model, i.e., a simula- 
tion without AGN heating. In this simulation we keep 
the central gas sink term but neglect feedback. This 
case leads to the classic 'cooling flow' catastrophe. The 
plasma in the central region of the cluster (~100-200 kpc) 
looses energy due to radiative cooling. Consequently, the 
gas pressure drops, leading to a subsonic spherical inflow. 
The increase in the central gas density results in higher 
cooling rates, which, in turn, accelerates the mass accre- 
tion rate, leading to runaway cooling and accretion. The 
results from a cooling flow run for the initial Tl-ratio of 
7 are presented in Figures [l] and [2j 

In the first column of Figure [l] we show the total mass, 
the average mass rate and the accretion rate as a func- 
tion of time. The quantities in the first two panels cor- 
respond to the gas within the central 20 kpc. Different 
gas phases are coded according to: 'cold' (T < 5 x 10 5 
K, blue), 'warm' (5 x 10 5 < T < 10 7 K, yellow), 'hot' 
(T > 10 7 K, red). The choice is not arbitrary because the 
limits of these temperature ranges correspond approxi- 
mately to different critical slopes in the cooling function 
A(T). Top panel reveals that the total mass of cold gas 
exceeds 10 12 M after a few Gyr. Middle panel shows 
that after 1 Gyr the cold phase continues to accumu- 
late in the core (20 kpc) at a steady average rate of 
~ 500 M Q yr" 1 (peaks can reach ~ 800 M Q yr" 1 ), a bla- 
tant discrepancy with observations. The cold gas never 
decreases with time, usually overwhelming the hot and 
warm phase by one and three orders of magnitudes, re- 
spectively. The accretion rate in the sink region is closely 
linked to the cooling rate. However, not all of the cooling- 
gas gets accreted because the inflowing cold gas forms a 
rotationally-supported structure. 

The second column displays the radial profiles of the 
emission-weighted temperature, electron number density, 
and entropy (for the X-ray emitting gas with T > 0.3 
keV). Curves are plotted every 500 Myr, and the color 
changes gradually from dark blue to light violet; black 
dashed lines correspond to the initial conditions. The 
pure cooling flow manifests itself as a sudden decrease in 
temperature, lesser than r v ; r /3. The projected emission- 
weighted (T ew ) temperatures falls in partic ular below 
the obser vational data o f Abell 1795 (circles , Tamura et 

1'he elec- 



al. 2001 



2002). 



and triangles, |Ettori et al. 
tron number density n e , linked to the surface brightness 
increases rapidly within 70 kpc from the center, with 
the simulated profiles lying above the observational data 
points. The entropy (fcbT/n 2 / 3 ) near 1 kpc is ~ 1 keV 
cm 2 , without showing any s ign of the typical observed 
floor riCavagnolo et al.||2008[ ). 

In Figure [51 we present the maps of X-ray surface 
brightness, the midplane cut of electron number density 
(n e ), and the gas phase diagrams with temperature ver- 




2 3 
time [Gyr] 



Fig. 1. — Diagnostics for the pure cooling flow simulation with 
initial minimum t coo i/tfl = 7. First column (from top panel): 
total mass and average mass cooling rate as a function of time in 
the central 20 kpc; accretion rate (in the sink region) as a function 
of time. Colors correspond to the three gas phases: blue - cold 
(T < 5 x 10 5 K), yellow - warm (5 X 10° < T < 10 7 K), red - 
hot (T > 10 7 K). Notice the strong cooling rate of ~ 500 M 
yr -1 . Second column: radial profiles of the emission-weighted 

temperature, electron number density, and entropy (k^T/n^ 3 ). 
Curves are plotted every 500 Myr, and the color changes gradually 
from dark blue to light violet; black dashed lines correspond to 
initial condit ions, which are base d on observat ional data of Abel l 
1795 (circles, [Tamura et al.|2001| and triangles, [Ettori et al.|2002| . 



sus n e (values are color-coded by gas mass in each bin). 
Each map has logarithmic scale and contains the data 
only for the central core (20 kpc); the snapshot is taken 
at 1.51 Gyr. The density maps reveal significant cold 
gas accumulation in the core, accompanied by a strong 
increase in the central surface brightness. Since in the 
'warm' temperature range the gas has a short cooling 
time, the ICM tends to accumulate in the 'cold' and 'hot' 
phases. This effect is clearly visible in the phase diagram. 

In the cooling flow run distinct cold gas clouds and fil- 
aments can be seen predominantly in the early stages of 
the evolution. At late times (e.g., 1.51 Gyr) the collapse 
proceeds in a monolithic way, and the atmosphere is glob- 
ally unstable. Note that the t coo i/tg criterion does not 
apply to a cooling flow, because the theoretical analysis 
assumes thermal balance (see Mil). In order to generate 
multiphase gas in a cooling flow, a substantial am plitude 
of perturbations is required, e.g., ^ 0.1 (see also Pizzo 



lato fc Soker||2005| Sll). Both initial Tl-ratio of 7 and 
21 (not shown) produce similar results. The latter case 
has longer saturation time and slightly smaller mass in 
the cold phase (the initial total mass in the box is also 
8% smaller). 

All the features of the pure CF model, discussed 
above, are in contrast with observations of galaxy 



clusters (Peterson et al. 2001 2003 Ettori et al. 2002 
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Fig. 2. — Pure CF simulation at 1.51 Gyr. Left: X-ray surface brightness (erg s _1 cm -2 ). Middle: n e midplane cut (cm -3 ). Right: 
joint distribution of logT (abscissa) versus logn e (ordinate), colored by gas mass in each bin (Mq); the size of the logarithmic bins is 
0.05. Each map has logarithmic scale (see color bars) and, apart from the joint distribution, is 20 kpc on a side (i.e., only the very central 
regions of the entire computational domain are shown). Thermal instabilities and extended cold gas are possible only in the early stage; 
later the collapse becomes monolithic, with a huge central accumulation of cold gas. 



Tamura et aL]|2003| |Peterson fc Fabian| |2006| . This 
simulation merely serves as a reterence case for the runs 
presented in the subsequent sections. A successful model 
will instead be characterized by a quasi equilibrium 
between heating and cooling, with low cooling rates 
(< 10% of the pure cooing flow), profiles consistent with 
observations (e.g., emission- weighted radial profiles with 
reasonable gradients and characteristic length scales). 
At the same time, the model needs to account for the 
spatially extended (up to ~ 15 — 20 kpc) distribution 
of cold blobs and filaments while not violating the 
constraints on the total mass in these features. Below 
we discuss models that meet these requirements. 



3.2. AGN Feedback: t coo i/tg = 7 

Before including AGN jet feedback, we verified that 
using the idealized prescription of Sll, i.e., local heat- 
ing rate given by the cooling rate averaged in spherical 
shells, produces thermal instabilities and spatially ex- 
tended multiphase gas. The results were consistent with 
Sll. However, unless cooling and heating are turned off 
in the central zones (as in Mil), this simple prescription 
led to the massive accumulation of cold gas in the cen- 
tral few kpc, especially when the gas dropout or outflow 
through a central boundary was neglected. These find- 
ings provide additional motivation to include a realistic 
feedback mechanism. Below we discuss the results of our 
simulations that include self-regulated feedback due to 
collimated AGN jets. 



3.2.1. 



6 x 10" 



Mechanical feedback due to AGN jets dramatically 
changes the entire evolution of the galaxy cluster, com- 
pared to the pure cooling case (Sec. 3.1). In Figures[3]and 
[4] we present diagnostic plots and maps for the runs that 
include AGN jets with mechanical efficiency e = 6 x 10~ 3 . 

The first column in Fig. [3] show the same quantities 
as those shown in the left panels in Fig. [TJ The gas 
mass in the cold phase M c (t), shown in the top panel, 



reaches ~ 10 11 M© at the end of the evolution. This 
is consistent with molecular cold gas observations (Edge 
2001 Salome et al.|2003 ), and is reduced by over a factor 
of 10 compared to the pure cooling flow value. Middle 
panel shows that the average cooling rate (M c ) is also 
significantly quenched: after an initial transition phase, 
in which cooling is slightly elevated, the system starts to 
self-regulate through more powerful outbursts, and even- 
tually reaches a steady rate of ~ 20 M yr _1 . This value 
is about 4% of the value seen in the pure CF simulation 
entirely cons istent with observations 



(.e.g. 



Peterson & 



Fabian||2006[ ). We note that the instantaneous cooling 
rate (not shown) is quite variable, changing rapidly from 
positive to negative values due to the action of AGN 
heating. 

As the cold gas mass M c in the central 20 kpc ex- 
ceeds the mass in the hot (3 times) and warm phases 
(10 3 times), the accretion rate is determined by the cold 
gas (last panel of first column); most of the time the 
black hole accretes cold gas at a rate between 1 — 5 M© 
yr -1 . Cold accretion rate (blue) exceeds hot accretion 
rate (red) by at least two orders of magnitude. The lat- 
ter is commonly associated with a Bondi-likc mode; since 
the hot accretion rate is so weak, feedback simulations 
based on Bondi prescription usually require an artificial 
boost factor or high mechanical efficienc ies > 0.1, even if 



the B ondi radius is resolved (e.g., Gila, Di Matteo et al. 



20081. However, we warn against interpreting our cold 



gas accretion rates literally, as the actual accretion rate 
onto the black hole will likely be lower. 

The mass accretion rate determines the AGN power 
and mass outflow rate (see Eq. ( 10 1 and (111). The prop- 



erties of one of the jets are shown in the second column 
in Figure |4j The cumulative energy injection as a func- 
tion of time is presented in the top panel. At 5 Gyr the 
total injected energy is ~ 10 62 erg. Even if the entire 
mass of the gas that falls into the sink region is added 
to the SMBH, its total mass will still be within reason- 
able limits (less than several 10 9 M©). Typical values 
of the jet power (second panel) and mass outflow rate 
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Fig. 3. — Diagnostics for the AGN feedback simulation with initial minimum t CO olAff = 7 and e = 6 X 10 — 3 (r7-6cm3). First column: 
cooling and accretion diagnostics; see the description of Figure^ for details. Second column (from top panel): total injected mechanical 
energy, outflow power and cumulative distribution function ('duty cycle') of the jet powers (100% corresponds to 5 Gyr); only for one of the 
jets. Third column: radial profiles; see the description of Figure^for details. Fourth column: radial profile of t coo \/tff and the temporal 
evolution of its minimum (gas with T > 0.3 kcV). Note that the cooling rate is reduced to ~ 4% of the pure CF value and that a steady 
average equilibrium is maintained for several Gyr. 
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Fig. 4. — Maps for the same simulation as in Fig. [3] (r7-6em3). See the description of Figure |2]for details. The maps correspond to the 
time of 1.51 Gyr since the beginning of the simulation. The salient feature is the amplification of thermal instabilities, due to the AGN 
perturbations, with cold blobs condensing out of the hot phase (t coo i/tff < 10) up to 10 kpc away from the center. 



are F jct ~ 3 - 4 x 10 44 erg s^ 1 and M jet - 0.5 - 1 M 
yr _1 . The cumulative distribution function (CDFpj; last 
panel) of the jet power suggests that the feedback is char- 
acterized by a 'duty cycle'. For example, assuming that 
the AGN is considered 'active' only when the mechan- 
ical power exceeds 10 44 erg s , the duty cycle would 



3 CDF(> X) represents the probability that the independent 
variable x is greater than value X; 100% value corresponds to 5 
Gyr (ttot). 



be about 70% (i.e., AGN active 70% of the time). This 
quantity may be poorly defined because of its sensitivity 
to the chosen threshold. In particular, deeper observa- 
tions may reveal more signs of AGN activity, leading to 
higher duty cycle s. Current observational estimates (e.g. 



Punnet al.|2006 ) place AGN duty cycle at a level > 70%. 

The radial profiles (Figure[3j third column) of the emis- 
sion weighted temperature, electron number density, and 
entropy reveal the presence of the feedback heating in 
two ways. First, the density inside 70 kpc is similar to 
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the initial (observed) conditions, meaning that the inflow 
due to the cooling is substantially quenched. Second, the 
temperature and entropy profiles do not fall below T v j r /3 
and 10 keV cm 2 , respectively. The profiles maintain a 
positive gradient, a sign that the cool core has not been 
destroyed by the outbursts. Notice that the real compar- 
ison with X-ray observations must be based on emission- 
weighted temperatures: the pure mass-weighted T profile 
would appear much more disturbed with higher central 
values. The electron number density profile (n e ), also 
linked to the X-ray emissivity, appears consistent with 
observational data, with a flatter gradient compared to 
the pure CF run. 

In Fig. (third column) we show the radial profile of 
icooiAff and its minimum as a function of time. Although 
we discuss the relationship between t coo \/ts and the pres- 
ence of the multiphase gas in more detail in Section 4, 
here we briefly point out that whenever this timescale 
ratio falls below pa 10, thermal instabilities grow and 
cold gas clumps condense out of the hot phase. Initially, 
the cold phase has a filamentary structure; then it turns 
into small clumps which fall toward the center in about 
a free-fall time. 

The SBx map and the the two-dimensional snapshot 
of n c at 1.5 Gyr, shown in Figure [4] (first and second 
panel), corresponds to a moment when the Tl-ratio falls 
below 10. The jet outburst creates two cavities along the 
z— direction, while at the same time the perturbations 
due to the jets seem to lead to the formation of dense cold 
blobs up to ~ 10 kpc away from the center. The presence 
of multiphase gas is also apparent in the n c — T phase 
diagram (Fig. Wl third panel) , which displays two distinct 
diagonal islands occupied by the gas in the cluster core, 
with small variations away from pressure equilibrium. 

Because the AGN heating is anisotropic close to the 
base of the jet and the cooling gas possesses small amount 
of angular momentum, the cold gas tends to accrete and 
condense more along the equatorial region. In a real- 
istic situation this gas cannot be entirely removed in a 
brief time (as in Sll); consequently, a cold rotationally- 
supported torus is formed, which provides most of the 
fuel for the AGN. The fuel supply from the torus, or 
the torus itself, can be temporarily disturbed by AGN 
outbursts, shocks, and the ICM turbulence. The cen- 
tral accumulation of cold gas could be further reduced 
by including star formation. 

Overall, the heating due to realistic AGN jets, with me- 
chanical efficiency ps 6 x 10~ 3 , appears to reproduce the 
main properties of the observed clusters: quenching the 
cooling flows while preserving the cool cores for several 
Gyr. At the same time, the model predicts the formation 
of the multiphase medium in the form of cold blobs and 
filaments while satisfying the o bservational constraints 
on the a mount of cold gas 
2011a|bp . 



(see 



McDonald et al. 2010 



3.2.2. e = 10 _i 

We now discuss a run for higher mechanical efficiency 
(e = 10 -2 ) compared to the one considered in the previ- 
ous Section (e = 6 x 10~ 3 ). We note that these values 
likely bracket the range of acceptable efficiency levels, 
as pointed out by previous computations (Gila; Sll) 
and observations (Merloni & Heinz 20081. We use the 



same set of diagnostic plots as before and we only dis- 
cuss the differences between these models. The results 
for e = 10~ 2 are shown in Figures [5l and |6j which are 
direct counterparts to Figures [3] and]4j respectively. 

The increase in the mechanical efficiency of the jets 
leads to stronger outflows characterized by a slightly 
more peaked jet power distribution centered around ~ 
6 x 10 44 erg s _1 . Because in the adopted feedback scheme 
Mj e t is proportional to the jet power, the outflows also 
reach slightly higher mass loading rates (> 0.5 M yr _1 ). 
This means that some jet events carry larger amount of 
linear momentum and, thus, have more 'piercing power' 
with which they plow through the ICM. Consequently, 
the jets release more of their energy at relatively larger 
distances from the cluster center, allowing a bit more 
cold gas accumulation along the equatorial region. The 
projected surface brightness map clearly shows the for- 
mation of a torus-like structure that was less pronounced 
in the lower efficiency run. 

The initial increase in the average cooling (and accre- 
tion) rate is strongly quenched at around 1 Gyr since the 
beginning of the simulation, which is earlier than in the 
previous run. After this short phase, the system gradu- 
ally settles to a quasi equilibrium state, with cooling rates 
around 40 — 50 M yr _1 , twice the level seen in the lower 
efficiency case. Paradoxically, more powerful jets are less 
efficient in stopping the cooling flow, due to tunnel carv- 
ing. The total injected energy is comparable in both 
runs (~ 10 62 erg), a sign that the self- regulation process 
is overall very similar, although with a slightly shorter 
active phase and a lower frequency of the outbursts (see 
jet power CDF). The total cold gas mass (2 x 10 11 M ) 
is slightly larger than in the preceding case, but still only 
10% of the pure cooling flow value, which is consistent 
with the observational constraints. 

During the entire evolution of this simulation, the out- 
flows tend to produce bigger bubbles or a deeper channel 
(see also the n e cut). The effect of more powerful jets 
is also to increase the scatter in the phase diagram: the 
gas in hot phase occupies a wider range of values in the 
n e — T plane about the (imaginary) diagonal line repre- 
senting average pressure. 

The minimum Tl-ratio falls below 10 only early in the 
evolution. Extended cold gas, in the form of blobs and fil- 
aments, is mainly formed at that epoch (see the n e cut). 
At later times, more powerful events raise i C ooiAff well 
over 10, and the ICM becomes hotter. Even if the system 
settles later around a minimum Tl-ratio of 13, extended 
cold filaments may be occasionally seen at moderate radii 
due to dredge up of the cold torus by outflows. The up- 
lift of the cold gas and the expansion of bubbles may 
further increase density, possibly triggering new thermal 
instabilities. This picture is also consistent with the ob- 
servational analysis presented in Mil (see their Fig. 12, 
right panel), suggesting that value 10 marks an approxi- 
mate transition to the regime where the multiphase gas 
may appear. 



3.2.3. 



6 x 10~ 3 , Wobbling 



To test the robustness of our results we included jet 
wobbling in the run with mechanical efficiency e = 
6 x 10~ 3 . Observational data is consistent with the pos- 
sibility that the jets may precess or their orientation may 
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Fig. 5. — Same as Figure [3] but for AGN feedback efficiency e = 10 -2 (r7-lem2). Raising the efficiency produces stronger and more 
massive outfiows, which have more difficulty in stopping the cooling flow, because they tend to release their energy at larger radii. Values 
are still consistent with observations. 
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Fig. 6. — Same as 
concentrated cooling. 



Figure U\ 
in the for 



but for AGN feedback efficiency e = 10 2 (r7-lem2). The time of the maps is 0.61 Gyr. Notice that 
5rm of a rotating cold torus, is favored over extended cold gas (i coo i/iff is usually > 10). 



change. Since the data is very limited, we keep the model 
very simple and change jet orientation whenever the ac- 
creted mass reaches a threshold of 10 7 M . For typical 
mass accretion rate of order ~ 1 M yr _1 , this corre- 
sponds to a reori entation timescale of ~ 10 Myr (e.g., 
Dunn et al.|2056| ). 

in the first simulation with jet wobbling, we assumed 
that the probability of jet pointing in a given direction is 
given by a random uniform distribution over An steradi- 
ans (the bipolar jets are always collinear). Initially, the 
jet orientation changes at a slower rate, ~ 100 Myr, be- 
cause the accretion rate is relatively low. Thus, in this 
early stage, the evolution is very similar to that for the 



non- wobbling case with efficiency of e = 6 x 10~ 3 . How- 
ever, the energy is better spatially distributed, and thus 
the cooling rate is more quenched, reaching values < 10 
M Q yr _1 . As the accretion rate accelerates, and the jet 
orientation starts to change on timescales of 10 Myr or 
less, the deposition of energy becomes nearly isotropic. 
A more uniform distribution of energy makes it more dif- 
ficult to feed the black hole via the cold torus, which is 
now strongly disturbed. Rather than being used to drill 
a channel through the ICM, the AGN jet is no longer able 
to efficiently deliver the energy to the ICM beyond r > 15 
kpc. This enables the gas to cool rapidly at these radii, 
which leads to the cooling catastrophe in about 1 Gyr. 
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Larger efficiencies could possibly stop the cooling, but 
the nearly isotropic centrally-concentrated energy distri- 
bution would destroy the cool core. 

In order to remedy this situation, we considered mod- 
els where the jet orientation is confined to two conical 
regions of half opening angle 75°, 45°, and 25°. As be- 
fore, the probability distribution of the jet orientations is 
randomly uniform within these regions. In the first two 
cases, the final result is similar to the case considered 
above, but the cooling flow is quenched for an additional 
time of < 1 Gyr. 

For half opening angle of 25°, the jet energy distri- 
bution becomes narrower and more anisotropic. Apart 
from larger bubbles and a more frequent jet channel frag- 
mentation, this computation resembles the previous best 
model r7-6em3: the cooling rate is steady for several Gyr 
(~ 10 — 30 M Q yr _1 ), with total cold gas mass around 
10 11 M.q. Cold blobs are produced mainly in the first 1.5 
Gyr when Tl-ratio oscillates around 10. At later times, 
the system settles to a hotter state characterized by the 
Tl-ratio of ~15-16. Jet events are also more variable 
than in the r7-6em3 model, changing rapidly from 10 43 
to 10 45 erg s _1 , because more inclined outbursts are able 
to partially dismantle the central torus, temporarily in- 
hibiting accretion. 

In summary, randomly wobbling jets with small and 
moderate opening angles produce similar results to 
unidirectional narrow outflows, especially for stopping 
the cooling flow catastrophe and generating extended 
multiphase gas. The other relevant feature is that 
the central accumulated cold torus can be frequently 
destroyed, depending on the wobbling parameters. 



3.3. AGN Feedback: t cool /t s = 21 

3.3.1. e = 6 x 1CT 3 

In order to study the effect of feedback in slower cooling 
systems, we generated new initial hydrostatic conditions 
as follows. We decreased the central density preserving 
the initial profile slope at r > 80 kpc; simultaneously, 
we increased the central temperature to ensure that the 
ICM remains in hydrostatic equilibrium. The new ini- 
tial conditions correspond to the minimum Tl-ratio of 
~ 21. The new initial profiles are plotted in Fig.[7j black 
dashed lines (third column). The entropy is elevated to 
30 keV cm 2 at 10 kpc - s till within observational bounds 



(Cavagnolo et al. 20091) - while the total gas mass is 

1 J -„ l — L-J tw — frn__ u. r. iL. 



decreased only by 87o. The results from the runs corre- 
sponding to these new initial conditions and AGN me- 
chanical efficiency ofe = 6xl0~ 3 are shown in Figures 
[7] and [8] These figures are analogous to those previously 
discussed. 

Overall, the behavior of the system is quite similar to 
what we observed before. The model is successful in 
predicting a globally stable ICM core, low cooling rates 
(~ 50 M Q yr _1 , 10% of pure CF), preservation of the 
cool core (e.g., entropy, SB X ), typical features of AGN 
jets (bubbles, shocks, turbulence), and the amount of 
multiphase medium consistent with observations. The 
main result is that the local thermal instability is trig- 
gered despite the fact that the initial minimum Tl-ratio 
exceeds ~ 10. It is only when the instantaneous value of 
this ratio falls below the critical threshold of ~ 10 that 



the multiphase medium suddenly appears in the simula- 
tion (e.g., after 1 Gyr). As before, this leads to increased 
feedback that brings the system back to the high Tl-ratio 
regime, which in this case is well above ~ 20 and higher 
than in the previously considered systems. 

Even if the typical jet power (few 10 44 erg s^ 1 ) and 
total injected energy (4 x 10 61 erg; single jet) are lower 
compared to r7-6em3, the AGN outflows have a relatively 
stronger impact on the lighter ICM, secularly increasing 
the Tl-ratio back above 20. As a consequence, spatially 
concentrated cooling, in the form of a rotating torus, 
is present during the entire evolution. In Section 4 wc 
make an in-depth comparison between the models with 
different initial t coo i/tg ratio. 



3.3.2. 



10" 



The total cold gas mass in the model discussed above 
reaches 2 x 10 11 Mq, which is higher than r7-6em3. Al- 
though this is observationally acceptable, some systems 
show smaller amounts of cold gas mass. Therefore, we 
experimented with increasing the mechanical efficiency 
in order to decrease the amount of cold gas (figures not 
shown). However, as argued in Sec. 3.2.2, higher effi- 
ciencies lead to stronger and more mass-loaded jets that 
deposit their mechanical energy farther away from the 
center. If in the initial outbursts the system appears 
to be more heated due to shocks, the situation reverses 
later, when the jet carves a deeper channel in the ICM. 
Thus, more accretion occurs along the equator and a 
strong cooling flow develops. As seen before, this feature 
could be alleviated by a randomly wobbling outflow, with 
a less anisotropic energy deposition. 

Overall, the simulations with higher minimum initial 
Tl-ratio 21, can again produce extended multiphase gas 
whenever the instantaneous t coo \/t$ falls < 10. In fact, 
in the early phase of the evolution, when cooling slightly 
prevails over heating, the Tl-ratio slowly decreases, cross- 
ing the threshold in about a Gyr. Later, stronger out- 
bursts set the system back to a hotter state and higher 
Tl-ratio, until the cycle restarts again (it may require 
several Gyr). The main trend for the entire 5 Gyr is 
therefore concentrated cooling occurring in the core. 

Contrary to the more massive models (Tl-ratio 7), 
the right efficiencies in order to quench the cooling 
flow and preserve the cool core should not be more 
than 6 x 10 -3 . A less massive system is in fact more 
susceptible to mechanical outbursts, and requires a 
more delicate feedback, in analogy to the environment 
of galaxy groups or ellipticals, as has been shown in 
previous works (Glib, Sll). 



4. DISCUSSION 

Galaxy clusters are extremely important astrophysical 
objects. Lying on the exponential tail of the halo mass 
distribution, they provide a sensitive probe of cosmol- 
ogy. They are also home to the most massive galaxies 
and black holes in the universe. Galaxy clusters are ob- 
served in all wavebands, from radio to X-ray and gamma 
rays. Owing to their deep potential wells, they are al- 
most closed systems with the baryon fraction similar to 
the universal value. Moreover, most of the baryons in 
clusters are in the form of hot X-ray emitting plasma, 
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Fig. 7. — Same as Figureplbut for initial minimum t coo i/% = 21 and AGN feedback e = 6 X 10 — 3 (r21-6em3). Starting with a higher 
Tl-ratio produces an initial phase of slow cooling, followed again by strong heating when the t coo i/% ratio approaches 10 (the multiphase 
gas threshold). AGN outflows have a relatively stronger impact on the lighter ICM, secularly increasing the Tl-ratio back above 20. This 
computation shows that the average thermal equilibrium is actually quasi stable, with phases of slight cooling or heating dominating 
throughout the entire evolution. 
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Fig. 8. — Same as [4] but for initial minimum i C ool/tff = 21 and AGN feedback e = 6 X 10~ 3 (r21-6em3). The time of the maps is 1.27 
Gyr. Notice that concentrated cooling is favored through the whole evolution, while cold extended blobs condense out of the hot phase 
only in short windows (again whenever t coo i/iff < 10). 



the intracluster medium. Understanding the state of the 
ICM will not only make clusters more useful cosmologi- 
cal probes, but also shed light on galaxy formation and 
on the state of most of the baryons in the universe. 

In this paper we try to understand the thermal state 
of the ICM core, closely tying our models with obser- 
vations. From multiwavelength observations it is clear 
that the dense plasma with short cooling times in cluster 
cores is not cooling in a catastrophic fashion. What is the 
mechanism at the base of such self-preservation behav- 
ior? Radio bubbles and X-ray cavities associated with 



cool-core clusters suggest that the active galactic nucleus 
acts as a thermostat, fueled by the cooling gas. As argued 
by Gila, the best realistic model able to maintain quasi 
thermal equilibrium for several Gyr, i.e., quenching the 
cooling flow without destroying the cool-core structure, is 
an asymmetrical input of mechanical energy, in the form 
of massive outflows. The cooling flow solution (Sec. 3.1) 
is very different from the observationally supported pic- 
ture where the cluster core is sustained by AGN heating. 
Mass accretion and cooling rates in our best feedback 
simulations (Sec. 3.2.1 and 3.2.2) are suppressed by a 
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factor > 10 compar ed to a pure cooling flow, in l ine with 



observation s (e.g., Peterson et al. 2001 2003 Tamura 



et al ||2003| |Petersonl?Tf bian 2006}. TEe X-ray sur 



face brightness is also not strongly peaked in the center; 
the entropy and emission-weighted temperatures do not 
rapidly decline under the initial values, conserving the 
positive gradient. Using mechanical efficiencies in the 
range 5 x 10 -3 — 10~ 2 , the typical outflow powers are 



~ 10 44 — 10 erg s , comparable to the X-ray luminos- 
ity. 

The fundamental result, consistent with observations, 
is that extended cold filaments and blobs condense out 
of the hot phase whenever t coo \/te ^ 10 or, equivalently, 
when the central entropy becomes < 30 keV cm 2 . The 
cold clumps fall then toward the center in about a free- 
fall time, forming a central cold rotating torus and in 
part boosting the black hole accretion rate. The for- 
mation of cold filaments is a result of interplay between 
thermal instability and gravity. Previous idealized anal- 
yses by Mil and Sll enforced strict thermal equilibrium 
and found that multiphase gas is generated according to 
the same thresholcjf] 

A realistic feedback displays a variegated self-regulated 
evolution, producing complex observed phenomena, 
which are not captured by simple idealized prescriptions. 
In the present work, we tested bipolar AGN jets, tied to 
the mass accretion rate at very small radiP] The asym- 
metrical input of energy produces strong elliptical shocks 
and regularly inflates hot bubbles, which rise outwards. 
At the same time the inflow is not entirely stopped: along 
the equatorial direction the gas, mainly cold, continues to 
fuel the black hole. Turbulence, entrainment and mixing 
dissipate the directional kinetic energy of the jets, further 
inhibiting the cooling flow near the center. When the jet 
ram pressure is substantial, the cold gas is also dredged 
up, as observed in a few cases such a s Perseus and Hydra 
(Hatch et al.||2006||Gitti et al.|2011b. The cycle repeats, 



when cooling starts to dominate again and the cold gas 
condensation ignites more powerful outbursts, setting the 
system in a slightly hotter state. 

Overall, thermal equilibrium is only roughly main- 
tained, on large spatial and temporal scales. For exam- 
ple, the two maps in Figure [9j for best model r7-6em3, 
display the net thermal energy increase (red) or decre- 
ment (blue) per volume and time. The internal energy 
equation is given by 



de 
M 



= -V • (ev) - P(V • v) - n n;A(T), 



(16) 



where e is the internal energy density (erg cm -3 ). The 
maps indicate that 'heating' (de/dt > 0) and 'cooling' 
(de/dt < 0) alternate on scales of a few kpc. Averages 
of de/dt in radial shells of width > 4 — 5 kpc suggest 
that a rough thermal equilibrium is maintained, although 
cooling appears slightly dominant in the central 10 kpc, 
while the jets tend to deposit their energy at larger radii. 
Equation ( 16 1 does not, however, fully capture the irre- 



4 The criterion in spherical coordinates is less stringent compared 
to the Cartesian case (~ 1) because instabilities are amplified by 
geometrical compression, as blobs sink toward the center. 

5 Our simulations do not resolve the sphere of influence of the 
SMBH, so the actual mass accretion rate might be smaller. Thus, 
in reality the feedback efficiency may be even larger. More work is 
needed to satisfactorily study the fate of the accreting cold gas. 



versible heating associated with turbulent mixing, shocks 
and, wave damping due to numerical effects. A better 
direct indicator of the balance between cooling (£) and 
heating (H) is the variation of the entropy S (oc lni^, 

where K = k^T/nl^ 3 ), which is insensitive to adiabatic 
compressions or expansions: 



PT 



dS 
~dt 



U-C, 



(17) 



where d/dt is the Lagrangian derivative. 

The right panels in Figure 9 show the time variations 
in radial shells of the 'astrophysical' entropy K for the 
whole ICM, with respect to the initial value Kq. As 
expected the entropy increases on a secular timescale 
due the continuous action of the AGN feedback. On 
smaller temporal scales, e.g., a few typical cooling times 
(100 — 200 Myr), the amplitude of the fluctuations in 
the six thick spherical shells (centered at r = 22, 18, 14, 
10, 6 and 2 kpc; shown from the top to bottom row, 
respectively) can reach 2 or 3 times the mean value, es- 
pecially in the nuclear region. Although the fluctuations 
in K can be substantial, its mean saturates at a stable 
level. This suggests that the cluster core does not dra- 
matically deviate from the state of (quasi) equilibrium 
between cooling and heating, at least on scales of tens 
kpc and hundreds of Myr. This was also suggested by 
the low cooling/accretion rates compared to a pure cool- 
ing flow. The assumption of strict global thermal equi- 
librium made in the work of Mil and Sll is therefore 
very idealized. It is thus particularly remarkable that 
our criterion for the formation of the multiphase gas is 
in good agreement with their result. 

Unlike in Sll, the central kpc region is not excised, 
allowing the cold gas (either nuclear cooling or fallen 
blobs) to form an angular momentu m-supported toru s 
as observed in the core of M87 (e.g., Ford et al. 1994). 
Such a disk can be partially dismantled when the jet di- 
rection is perturbed by turbulence or by random jet pre- 
cession (Sec. 3.2.3). Nevertheless, the central dens e disk 
is a common feature among elliptical galaxies (e.g., Mac 



I chetto et aLpJ97] [Mathews fc Brighenti|2003|| Verdoes et 



|al.||2006p . Central y concentrated atom ic gas i s also seen 
in the survey of|McDonald et al^pOToj [20rTa| . The cen 



tral disk is certainly a nursery tor new stars. Since the 
star formation efficiency is quite uncertain (few per cent, 
up to 50%, with an av erage value of 15%, according to 



McDonald et al. 2011b I we decided to not model star for- 
mation via a sink term. Thus, the central accumulation 
of cold gas could be further reduced. 

The presence of multiphase gas is clearly seen in the 
n c —T joint distributions. The cold dense phase (10 4 — 10 5 
K), in the upper left corner, is present in large quantities 
and detached from the hot rarefied gas (> few 10 7 K), 
with a lack of g as at intermediate temperatu res. A re- 
cent Ha survey ( McDonald et al.|2010 2011a) has shown 
that extended cold filaments are found in > 35% of all 
clusters. It is striking that some of their maps (see their 
Fig. 4) resemble our computed morphologies, with ex- 
tended blobs and filaments up to ~ 20 kpc: e.g., Abell 
0496, Hydra A, Abell 1644, Sersic 159-03, and Abell 1795 
(the cluster on which our model is based). Another case 
with substantial extended Ha emission is the Perseus 



cluster (Salome et al. 2006 Fabian et al. 20081. More 
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x [kpc] time [Gyr] 

Fig. 9. — Thermal equilibrium diagnostics (best model r7-6em3). Left: 2D midplane cut of the cluster core (20 kpc) showing the net 
internal energy increase (top) or decrement (bottom), averaged between 1 and 2 Gyr (erg s _1 cm -3 ; in logarithmic scale). Right: temporal 
variation of entropy of the whole ICM with respect to the initial value Kq. The six panels show (K — Ko)/Ko as a function of time, every 
10 Myr, at different spherical shells of width 4 kpc centered at r = 22, 18, 14, 10, 6 and 2 kpc (from top row). 



over, the total cold gas mass we obtain, ~ 1 M , is from blue (early) to red (late). For the run with the 

initial t coo i/tg = 7 (upper panel), the cold gas is much 



consistent with the constraints of Edge (2001 ), whose ob- 
servations indicate clusters harboring up to 3 x 10 11 M 
of molecular gas. 

The white noise perturbations in the initial conditions 
are wiped out over a cooling time. In fact, in the pure 
cooling flow computation, multiphase gas is only possi- 
ble in such restricted interval, while the secular evolu- 
tion is dominated by monolithic cooling. On the con- 
trary, AGN feedback drives and sustains the fluctuations 
in the long term evolution, although the original seeds 
initially help to distribute better the jet energy due to 
the interaction of jet with the inhomogeneous ICM. The 
thermal state of the cluster core depends mainly on the 
feedback efficiency: with smaller e the frequency of cold 
filaments condensing out of the hot phase is larger (Figs. 
3, 5, 7). Nevertheless, higher initial t coo \/tg (and lower 
density) generates extended multiphase gas according to 
the TI criterion, but thermal equilibrium is reached at 
later epochs. It is interesting to note that AGN feedback 
is more effective in lighter systems, requiring therefore 
lower efficiencies (< 6 x 10~ 3 ), as has been shown for 
galaxy groups (Glib, Sll). 

To better quantify the presence of extended multiphase 
gas, Figure[T0]shows a scatter plot of the average distance 
of cold gas (-RcoidS excluding r < 3.5 kpc) as a function of 
tcooi/ts- The points are color coded according to time: 



more spatially extended, especially in the first 2 Gyr. 
The run with less dense initial conditions (Tl-ratio 21) 
shows a lower frequency of cold blobs at large distances 
because t coo \/ts > 10 most of the time (compare Figs. 3 
and 7). In both cases, the points clustered around ~ 4 
kpc correspond to the rotating cold torus. Since this plot 
shows the current average distance rather than the radius 
where the cold gas originates, extended filaments appear 
even when Tl-ratio is elevated to ~ 12 — 15 due to AGN 
heating. This tendency is further enhanced if cold blobs 
are dredged up by powerful outflowing jets. Moreover, 
as the Tl-ratio is averaged in spherical shells, a few cold 
zones may contribute to i? C oid ; without changing much 
the average t coo \/tg. That is, some extended Ha emission 
is consistent with t coo \/tg slightly greater than 10. 

Figure 11 displays i? C oid gas as a function of time. With 
denser initial conditions (r7-6em3), i C ooiAff is less than 
10 for several short intervals within 2 Gyr and we expect 
extended cold ga^] Indeed, the average distance of the 
cold gas is significantly larger than that of the cold torus 
(« 4 kpc) during this period. Similarly, for the lower ini- 
tial density run we see extended cold filaments at 1 — 1.5 
Gyr, again when t coo \/tf[ < 10; later times are however 

6 Note also that Af acc is boosted exactly during this period (see 
Fig.[3f. 



14 



Gaspari, Ruszkowski & Sharma 



12 



en 6 







5 10 15 


20 




m in (tccoiAf) 














6 




Fig. 10. — Extended cold radius as a function of minimum 
shell-averaged Tl-ratio. This radius is defined as the average 
distance from the center of cold cells that do no belong to the 
bulk of concentrated cooling (r < 3.5 kpc). The points are 
color coded according to time (Gyr; see color bar) Up: r7-6em3 
simulation. Down: r21-6em3 simulation. Black line represents 
the linear regression fit to the scatter plot (points are drawn 
every 10 Myr). Notice that these plots only show where the cold 
gas instantaneously resides (on average), not where it originally 
formed. 



dominated by the nuclear cold gas. We conclude that 
the instantaneous value of t coo i/tg is a robust indicator 
of extended cold gas formation, while nuclear cold gas 
can be prese nt i rrespective of the Tl-ratio (e.g., bottom 
panel of Fig. 10). 

It is likely that other physical processes play an im- 
portant role in shaping the morphol ogy of condensing 



multiphase gas. Sharma et al. (2010) and Mil showed 
that anisotropic thermal conduction leads to cold gas 
stretched along the local magnetic field direction. In 
order to be able to fully understand the exact topol- 
ogy of cold filaments, this physics should be included 
in the simulations (with an adequately higher resolu- 
tion). However, the generation of cold filaments, when- 
ever t coo \/tg < 10, is not likely to be affected by the 
presence of anisotropic thermal conduction. 

The other physics that is missing in our simulations 
involves cosmological accretion and mergers, which drive 
additional turbulence, inducing more instabilities. Like 
AGN feedback, major mergers are associated with en- 
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Fig. 11. — Extended cold radius as a function of time. Up: 
r7-6em3. Down: r21-6em3. Starting with lower Tl-ratio leads to 
an evolution characterized by both extended (t coo i/% < 10) and 
concentrated cold gas. On the contrary, extended multiphase gas 
is rarer in the opposite starting regime, although not absent. 



tropy generation and strong heating of cluster cores, rais- 
ing the cooling time over the cluster age. Thus, mergers 
may be responsible for extreme non-cool-core clusters. 
However, they are also rare with timesca les much longer 



than the cluster core cooling times (e.g., Fakhouri et al. 
20101. Therefore, AGN feedback, rather than mergers, is 



expected to be responsible for balancing cooling in most 
systems. 

Our theoretical models are able to solve the cooling 
flow problem for several Gyr. The key factor resides 
in adopting realistic mechanical AGN jets. The bulk of 
AGN heating is in fact a combination of shock dissipation 
and g radual thermalizatio n of the mechanical energy (see 
also Briiggen et al. 2007). Random motions, either due 
to AGIN, multiphase gas, and cosmological turbulence, 
strongly help the deposition of the feedback energy in 
the cool core , in contrast to idealized sy mmetric atmo- 
spheres (e.g., Vernaleo k. Reynolds||2006 1 . 

An important issue is whether our results depend 
sensitively on resolution. The morphology of the cold 
phase depends on the resolution. We note that finite 
resolution may prevent the disruption and mixing of 
the small infalling clouds with the ambient ICM. Thus, 
limited resolution could in fact mimic some of the effects 
of the physics that we do not include (e.g., magnetic 
draping). We performed a test with two times lower 
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resolution, but with the same resolution in the jet 
injection region (r7-6em3c), and obtained very similar 
results. In particular, a steady global thermal equilib- 
rium is reached in a few Gyr, with only a slightly larger 
final cold mass (2 x 10 11 M Q ). The evolution of the 
Tl-ratio and condensed cold blobs appears also similar 
to the higher resolution model (r7-6em3), although the 
increased numerical dissipation tend to slow down the 
growth of thermal instabilities at larger radii. We note 
also that our simulations have a resolution ~ 10 times 
higher compared to the work done in Gila. Although 
in Gila we adopted a cold mass dropout and a different 
jet injection method, AGN jets could also stop the 
cooling flow, in line with observations. Long term 
thermalization seems therefore not dramatically affected 
by resolution or the exact heating prescription. The fine 
details of the feedback can however change with higher 
resolution, e.g., jets tend to become more frequent (cf. 
Gila). 



5. CONCLUSIONS 
Our main conclusions can be summarized as follows: 



Mechanical AGN feedback self-consistently linked 
to BH accretion (with efficiency 5 x 10~ 3 — 10 -2 ), 
can solve the cooling flow problem, quenching the 
cooling rate down to ~ 10 M yr _1 (1-10% of 
the pure cooling flow value), preserving the cool 
core structure in agreemen t with observations 
dPeterson et al.| |2001[ j2003 El urn el 



i'amura et al.f 
McNamara & 



HI [2002 



Vulscn 



powers are 
~ 10 4 km s 



,44 



10 



J eterson &; Fabian 2006 



2007) 



erg s" 



Typical outflow 
_1 , with velocities 
; rates ~ M ffl yr -1 



(see 


Crenshaw et al.||2003[ |Morganti et al. 


2005 


2007 


N csvad ba et al.||2008l |20Tl| 


Tombesi et al. 



with fixed inclinations or fast random wobbling 
with half-opening angle < 25°. 



Multiphase gas, with cold (T ~ 10 4 K) and hot 
(T > 10 7 K) components, is a direct consequence 
and source of feedback. The cold gas comes in 
two phases: extended phase (up to 10-20 kpc of 
distance form the center) in the form of blobs and 
filaments generated by thermal instabilities (and 
occasionally due to dredge-up by jet ram pressure) 
and a more spatially concentrated one in the form 
of central rotating torus due to nuclear cooling 
or accumulated fallen blobs. The Ha survey of 

shows both kinds 



McDonald et al. (2010 2011a 



of Ha emitting gas. Nuclear cold reservoirs of gas, 
especially in the form of rotating dis ks are also 
a common feature among ellipticals (|Ford et al. 
1994} [Macchetto et al ||T997| patliews & BrignentT 
2003; Vcrdoc s et al.||2006| a star formation sink 
term and jet wobbling can reduce it considerably 
though). Total cold gas masses are ~ 10 11 M Q 
after 5 Gyr, within the observational constraints 
(< 3 x 10 11 M , |Edge||200l| |Salome et~al]|2003 1 . 



• The generation of extended cold gas is dictated 
by an important threshold: whenever the instan- 
taneous ratio of the cooling time to the free-fall 
time is < 10 thermal instabilities lead to cold gas 
condensing out of the hot ICM at large radii. On 
the contrary, when t coo \/tg > 10, slow runaway 
cooling occurs only within the central few kpc. 
This threshold seems to be surprisingly robust, as 
predicted by simple analytical arguments (Mil) 
and simulations with idealized feedback (Sll). The 
criterion is also supported by recent observations 
(Fig 12 in Mil; |Voit et al.||2008| |Cavagnolo et al. 



2008) 



• Multiphase ICM can exist in cool-core clusters irre- 
spectively of their initial minimum t coo i/ts. With 
a realistic feedback the cluster core always finds a 
way to self-regulate toward a state of global quasi 
thermal equilibrium on large spatial and temporal 
scales. The system forgets about the initial state 
in a few cooling times and the eventual evolution 
depends primarily on the feedback efficiency. 
Lower e results in frequent episodes of extended 
cold gas. Intense cooling episodes are followed by 
increased jet activity and the core returns to a 
higher entropy level. This phase is more prolonged 
and intense in less massive systems (e.g., initial 
icooiAff = 21). The reduced gas inflow leads to 
the decreasing jet power, allowing a slow spatially 
concentrated cooling to start again (t coo i/tg > 10). 

Overall, the main thermal properties of cluster cores 
are determined by global quasi-thermal equilibrium, 
with an amazingly detailed complexity of cold fila- 
ments/torus, jets, bubbles and shocks. Additional 
processes, such as thermal conduction, magnetic fields 
and cosmological accretion, are required to make more 
precise comparisons with observations. Nevertheless, 
we showed that multiphase gas due to thermal insta- 
bilities is an essential element of the feedback cycle: 
end product of cooling and fuel for the central black hole. 

ACKNOWLEDGMENTS 

The software used in this work was in part developed 
by the DOE NNSA-ASC OASCR Flash Center at the 
University of Chicago. We acknowledge the NASA 
awards SMD-10-1609, SMD-11-2209 (Pleiades), the 
CINECA awards HP10BPTM62, HP10BOB5U6 (SP6), 
and the University of Michigan for the availability of 
high performance computing resources and support. MG 
would like to thank F. Brighenti for useful discussions 
and P. Temi, as the principal support scientist at the 
NASA/ Ames base. MR acknowledges the NSF grant 
1008454. Support for PS was partially provided by 
NASA through Chandra Postdoctoral Fellowship grant 
PF8-90054 awarded by the Chandra X-Ray Center, 
which is operated by the Smithsonian Astrophysical 
Observatory for NASA under contract NAS8-03060. 
We are also thankful for the hospitality of the Kavli 
Institute for Theoretical Physics (KITP) at UC Santa 
Barbara, where this work was initiated. 



1G 



Gaspari, Ruszkowski & Sharma 



REFERENCES 



Balbus S.A., Soker N., 1989, ApJ, 341, 611 

Blanton E.L., Sarazin C.L., McNamara B.R., Wise M.W., 2001, 
ApJ, 585, 227 

Bogdanovic, T., Reynolds C.S., Balbus S.A., Parrish I.J., 2009, 
ApJ, 704, 211 

Bohriiiger H., Voges W., Fabian A.C., Edge A.C., Neumann 

D.M., 1993, MNRAS, 264, L25 
Brighenti F., Mathews W.G., 2002, ApJ, 567, 130 
Brighenti F., Mathews W.G., 2003, ApJ, 587, 580 
Briiggen M., Heinz S., Roediger E., Ruszkowski M., Simionescu 

A., 2007, MNRAS, 380, 67 
Cavagnolo K.W., Donahue M., Voit G.M., Sun M., ApJ, 683, 107 
Cavagnolo K.W., Kenneth W., Donahue M., Voit G.M., Sun M., 

ApJS, 182, 12 

Crenshaw D.M., Kraemer S.B., George I.M., 2003, ARAA, 41, 117 
Conroy C, Ostriker J.P., 2008, ApJ, 681, 151 
Di Matteo T., Colberg J., Springel V., Hernquist L., 2008, ApJ, 
676, 33 

Dunn R.J.H., Fabian A.C., 2006, MNRAS, 373, 959 
Edge A. C, 2001, MNRAS, 328, 762 

El-Zant A.A., Kim W.T., Kamionkowski M., 2004, MNRAS, 354, 
169 

Ettori S., Fabian A.C., Allen S.W., Johnstone R.M., 2002, 

MNRAS, 331, 635 
Fabian A.C., 1994, ARAA, 32, 277 

Fabian A.C., Sanders J.S., Crawford C.S., et al., 2003, MNRAS, 
344, L48 

Fabian A.C., Johnstone J.S., Sanders J.S., et al., 2008, Nature, 
454, 968 

Fakhouri O., Ma C, Boylan-Kolchin M., 2010, MNRAS, 406, 2267 
Finoguenov A., Jones C, 2001, ApJ, 557, L107 
Finogucnov A., Ruszkowski M., Jones C, ct al., 2008, ApJ, 686, 
911 

Ford, H. Harms, R.J., Tsvetanov Z.I., et al., ApJ, 1994, 435, 27 
Fryxell B., Olson K., Ricker P., et al., 2000, ApJS, 131, 273 
Gaspari M., Mclioli C, Brighenti F., D'Ercole A., 2009, in Heinz 

S., Wilcots E., eds., AIPC 1201, 'The Monster's Fiery Breath: 

Feedback in Galaxies, Groups, and Clusters', 309 
Gaspari M., Melioli C, Brighenti F., D'Ercole A., 2011a, 

MNRAS, 411, 349 [Gila] 
Gaspari M., Brighenti F., D'Ercole A., Melioli C, 2011b, 

MNRAS, 415, 1549 [Glib] 
Gitti M., Nulsen P.E.J, David L.P., McNamara B.R., Wise M.W., 

2011, ApJ, 732, 13 
Guo F., Oh S.P., 2008, MNRAS, 384, 251 
Guo F., Oh S.P., Ruszkowski, M., 2008, ApJ, 688, 859 
Hatch N.A., Crawford C.S., Johnstone R.M., Fabian A.C., 2006, 

MNRAS, 367, 433 
Jones C, Forman W., Vikhlinin A., et al., 2002, ApJ, 567, L115 
Macchetto F., Marconi A., Axon D.J., Capetti A., Sparks W., 

Crane P., 1997, ApJ, 489, 579 
Mathews W.G., Brighenti F., 2003, ARAA, 41, 191 
McCourt M., Sharma P., Quataert E., Parrish I., 2011, preprint 

(astro-ph/1105.2563) [Mil] 
McDonald M., Veilleux S., Rupke D.S.N., Mushotzky R., 2010, 

ApJ, 721, 1262 

McDonald M., Veilleux S., Mushotzky R., 2011a, ApJ, 731, 33 
McDonald M., Veilleux S., Rupke D.S.N., Mushotzky R., 
Reynolds C, 2011b, ApJ, 734, 95 



McNamara B.R., Nulsen P.E.J. , 2007, ARAA, 45, 117 
Merloni A., Heinz S., 2008, MNRAS, 388, 1011 
Morganti R., Tadhunter C.N., Oosterloo T., 2005, A&A, 444, L9 
Morganti R., Holt J., Saripalli L., Oosterloo T., Tadhunter C.N., 

2007, A&A, 476, 735 
Morsony B.J., Heinz S., Briiggen M., Ruszkowski M., 2010, 

MNRAS, 407, 1277 
Narayan R., Medvedev M.V., 2001, ApJ, 562, L129 
Navarro J.F., Frenk C.S., White S.D., 1996, ApJ, 462, 563 
Nesvadba N.P.H., Lehnert M.D., De Breuck C, et al., 2008, 

A&A, 491, 407 

Nesvadba N.P.H., De Breuck C, Lehnert M.D., et al., 2011, 
A&A, 525, 43 

O'Dea CP., Baum S.A., Privon G., et al., 2008, ApJ, 681, 1035 
Parrish I.J., Quataert E., Sharma P., 2010, ApJ, 712, L194 
Peterson J.R., Paerels F.B.S., Kaastra J.S., ct al., 2001, A&A, 
365, L104 

Peterson J.R., Kahn S.M., Paerels F.B.S., et al., 2003, ApJ, 590, 
207 

Peterson J.R., Fabian A.C., 2006, PhR, 427, 1 
Pizzolato F., Soker N., 2005, ApJ, 632, 821 

Raffcrty D.A., McNamara B.R., Nulsen P.E.J., Wise M.W., 2006, 
ApJ, 652, 216 

Raffcrty D.A., McNamara B.R., Nulsen P.E.J., 2008, MNRAS, 
687,899 

Ruszkowski M., Begelman M.C., 2002, ApJ, 581, 223 
Ruszkowski M., Briiggen M., Begelman M.C., 2004a, ApJ, 611, 
158 

Ruszkowski M., Briiggen M., Begelman M.C., 2004b, ApJ, 615, 
675 

Ruszkowski M., EnBlin T.A., Briiggen M., Heinz S., Pfrommer 

CM., 2007, MNRAS, 378, 662 
Ruszkowski M., EnBlin T.A., Briiggen M., Begelman M.C, 

Churazov E., 2008, MNRAS, 383, 1359 
Ruszkowski M., Oh S.P., 2010, ApJ, 713, 1332 
Ruszkowski M., Oh S.P., 2011, MNRAS, 414, 1493 
Salome P., Combes F., 2003, A&A, 412, 657 
Salome P., Combes F., Edge A.C, et al., 2006, A&A, 454, 437 
Scannapieco E., Briiggen M., 2008, ApJ, 686, 927 
Sharma P., Chandran B., Quataert E., Parrish I. 2009, ApJ, 699, 

348 

Sharma P., Parrish I., Quataert E., 2010, ApJ, 720, 652 
Sharma P., McCourt M. , Quataert E., Parrish I., 2011, preprint 

(astro-ph/ 1106.4816) [Sll] 
Sutherland R.S., Dopita M.A., 1993, ApJS, 88, 253 
Tamura T., Kaastra J.S., Peterson J.R., et al., 2001, A&A, 365, 

L87 

Tamura T., Kaastra J., Makishima K., Takahashi I. 2003, A&A, 
399, 497 

Tamura T., Kaastra J.S. , den Herder J.W.A., Bleeker J. A.M., 

Peterson J.R., 2004, A&A, 420, 135 
Tombesi F., Cappi M., Reeves J. N., et al., 2010, ApJ, 521, 57 
Townsend R. H. D., 2009, ApJS, 181, 391 

Verdoes K., Gijs A., van der Marel R.P., Noel-Storr J., 2006, AJ, 
131, 1961 

Vernaleo J.C., Reynolds C.S., 2006, ApJ, 645, 83 
Voit CM., Cavagnolo K.W., Donahue M., et al., 2008, ApJ, 681, 
L5 

ZuHone J.A., Markevitch M., Johnson R.E., 2010, ApJ, 717, 908 



